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^p^l Abstract. We study the surface plasmon modes of an arbitrarily shaped nanoparticle 

• in the electrostatic limit. We first deduce an eigenvalue equation for these modes, 

, expressed in terms of the Dirichlet-Neumann operators. We then use the properties of 
these pseudo-differential operators for deriving the limit of the high-order modes. 
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1. Introduction 

The interaction between light and metallic nanoparticles can become very strong due to 
the excitation of surface plasmons. These hybrid modes of the electromagnetic field and 
the electron gas are confined to the surface of the particle and give rise to an enhance- 
ment of the incident field by several orders of magnitude [H El [3] . This enhancement 
enables a variety of applications ranging from the well-established surface-enhanced Ra- 
man spectroscopy (SERS), which allows the detection of even a single molecule [HE], 
to the emerging field of plasmonics [6], [7] , including for instance plasmonic waveguides 
which effectuate optical energy transfer below the diffraction limit [HI [8], [9]. 

While an exact analytical description of the optical response of a metallic nanopar- 
ticle exists only for very specific geometries like a sphere or an ellipsoid, numerical 
methods can be applied for a particle with an arbitrary, realistic shape. There is a wide 
range of numerical methods for light scattering [10]; the finite difference time domain 
approach (FDTD) is a very common one (originally proposed by Yee [H]). FDTD com- 
bined with a suitable discretisation of the particle enables one to calculate the response 
of an almost arbitrary particle to nearly any incident field. Moreover, in the case of a 
sufficiently small particle, i.e. for particles that are much smaller than the significant 
wavelengths, for which the dielectric function can be taken as constant, the surface 
plasmon resonances can be determined in the electrostatic limit via the eigenvalues of a 
surface-integral operator [121 [131 [E] • Hohenester and Krenn [H] have shown that this 
boundary integral approach can be used to calculate the surface plasmon resonances of 
single and coupled spheres, cylinders and cubic-like nanoparticles. 

In this paper we consider an arbitrarily shaped nanoparticle (see Fig. [1]) and show 
that the surface plasmon resonances in the electrostatic limit may be obtained as eigen- 
values of an operator which may be represented in terms of Dirichlet-Neumann op- 
erators. We argue that this reformulation makes the well-developed analytical tools 
pertaining to the study of such operators available for the study of surface plasmon res- 
onances in cases where exact solution formulae are not available. To support this claim 
we show here how the asymptotic behaviour of the high-order surface resonances may be 
analysed. To this end we use the well-known fact that the Dirichlet-Neumann operators 
are pseudo-differential operators, and use standard properties of such operators. 

In Section [2] we reformulate the boundary value problem as an eigenvalue problem 
and in Section [3] we show how one can apply the formalism to a half-space. In Section H] 
we introduce pseudo-differential operators and prove in rigorous mathematical terms 
the convergence of the high order modes. 
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Figure 1. Sketch of an example geometry, with quantities inside the particle being 
labelled by the index - and outside by +. S denotes the surface of the particle. 

2. Eigenvalue equation 

In the electrostatic limit the surface plasmon resonances of a particle as sketched in 
Fig. [1] are characterised by nontrivial solutions of the Poisson equation without external 
charges, 

A0±(r) = for r ^ 5, (1) 

(/,_(r) = </)+(r) hires, (2) 

edn(j)-{r)=dnMr) for res, (3) 

with the potential 0_ and 0+ inside and outside the particle, respectively; e is the per- 
mittivity of the particle and dn is the outward normal derivative on the surface S of the 
particle. For convenience we assumed by stating the boundary condition in eq. ([3]) that 
the particle surrounded by vacuum is homogeneous, isotropic, local and non-magnetic. 
Since we restrict ourselves to surface modes we further assume that the potential van- 
ishes for r ^ oo, so that is determined uniquely. We assume that the surface S is 
smooth, i.e. has no singularities. So surface resonances are given by thus values of e for 
which nontrival solutions of the system ([I])-© exist. 

In order to recast the problem dH)-© as an eigenvalue problem we introduce 
the Dirichlet-Neumann operators {D_ and D^) inside and outside the particle. The 
Dirichlet- Neumann operator inside the particle is defined as 

D_ : C°^{S) C°^{S) 

f ^ dn(j)f, 

where 0/ is the solution of the Dirichlet problem, i.e. it satisfies A0j(r) = inside and 
4>f{r) = f{r) on the surface 5* of the particle. The Dirichlet-Neumann operator outside 
the particle (-D+) is defined in a similar way. 

With the help of the operators D_ and the problem in eqs. ([I])-© can be recast 
in a compact manner as 

eD_f = D^f (4) 



Surface plasmon resonances of an arbitrarily shaped nanoparticle 



4 



where / is the restriction of 0± to S. Nontrivial solutions of this equation can be found 
for such e for which 

ker(eL>_ - D+) ^ ^ ker(eL)_L);^ -1)^0. (5) 

Note that is an invertible operator. Thus, the desired values of e are the inverse 
eigenvalues of the operator D^D^^, i.e., 



eigenvalue {D_D_^_^} 

Hence, by means of the Dirichlet- Neumann operators we have reformulated the 
boundary value problem for the determination of the surface modes as an eigenvalue 
problem. 

As a side remark, note that D_D^^ is not selfadjoint with respect to the standard 
scalar product, although both D±. are. This can be remedied either by considering 
the non-standard scalar product {f,g) = fg{D^f){x)g{x) dS{x) instead or by using 
D^^^^D_D^^^^ in ([5]) instead of D_D^^, which is selfadjoint for the standard scalar 
product and therefore shows the reality and completeness of the spectrum. The 
discreteness of the spectrum of D^D^^ is less obvious and will be shown below for 
a bounded surface S. 



3. Surface modes of a half-space 

As an illustration of the method developed above we determine the surface modes of a 
half-space. Since the resonances can be found at such e which are equal to the inverse 
eigenvalues of D^D^^ we first derive the Dirichlet- Neumann operators D_ and for 
a half-space. The surface S of the material, located at 2 < with the permittivity e, is 
chosen to be the x-y-plane, so that we have to solve the Laplace equation 

A<p^{x,y,z) = (7) 

for z < and 2; > 0, resp., with the Dirichlet boundary condition 

(pTi^, y, 0) = f{x,y); (8) 

outside the material, for 2; > 0, we assume vacuum. In order to solve eqs. (I?!) and (IE]) 
we use the Fourier transform in the x- and y-coordinates, 

f{tz) = j<lxe-'^<f{x,z), (9) 

with the two Fourier variables ^ and x = {x, y) . Fourier transformation of eqs. ([7j) and 
([8]) gives 

+ dl^^it z) = and 0^(^,0) = /(|). (10) 
From this one obtains the solution 



(|,;.) = e±l«IV(0, (11) 
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(using the additional boundary condition that vanishes for \z\ oo) from which it 
becomes obvious that the surface modes are confined to a small neighbourhood of the 
surface of the material. 

For the case considered the outward directional derivative is simply given by dz, so 
that the Dirichlet-Neumann operator can be stated as 



In particular, for a half-space we have D_ = — -D+ and therefore D^D^ = —1. Since 
the only eigenvalue of the unit operator 1 is 1 we obtain the well-known result that 
e = — 1 for the surface modes of a half-space. 

4. High-order surface modes 

Before we study the convergence of the eigenvalues of the operator D^D^^ in rigoruos 
mathematical terms, we briefly discuss the physical expectation. For example, 
considering a sphere surrounded by vacuum it is well known that the resonances of 
the surface modes are given by [15\ 



for k = 1,2,.... Obviously, for = 1, as corresponding to the dipole mode of the 
sphere, the resonance can be found at e = —2. On the other hand, in the limiting case 
of — s> oo, corresponding to the multipole mode of infinite order, the resonance occurs 
at e = —1, i.e., one retrieves the half-space result deduced above. This convergence 
property can be understood from the physicist's point of view as follows: The multipole 
modes of order k ^ 1 are related to fields which vary on length scales much smaller than 
the radius of the sphere; the higher the order k the smaller the length scale. Therefore, 
these high-order modes cannot distinguish the sphere from a half-space, since the sphere 
can be considered as locally fiat on the length scale of these modes. Hence, the con- 
vergence to e = — 1 for — > oo is not restricted to spheres, but should be expected for 
arbitrary particles. In the following we prove this property in strict mathematical terms. 

To this end we first briefly introduce the concept of pseudo-differential operators, 
for a more detailed introduction to this topic we refer the reader to [16] . The definition 
of a general pseudo-differential operator P with symbolp{x, ^) operating on the function 
^(a:;) can be stated as (with x,^ E R") 




(12) 



and as 




(13) 




(15) 
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From this definition it can be seen that if the symbol p{x, ^) is a polynomial in ^ the 
operator P is a conventional differential operator, which has constant coefficients if p is 
independent of x. However, for pseudo-differential operators one admits more general 
symbols: The smooth function p is required to have an asymptotic expansion as ^ — > oo 
of the form p{x,$,) ~ pm{x,$,) +pm-i{x,$,) +pm-2{x,$,) H , where each p^-i is posi- 
tively homogeneous in £ of degree m — i, i.e. satisfies Pm-i{x,$,) = \^\''^~''Pm-i{x, ^) for 
^ 7^ 0, where ^ = The number m can be an arbitrary real number and is called 
the order of P. The leading term pm is called the principal symbol of P and denoted 
as am{P)- In the case of a differential operator with p{x,$^) = ^ a^(a;)^^ (where 

/3 G Nq is a multi-index), the principal symbol is crm{P) = 

|/3|=m 

We defined pseudo-differential operators as acting on functions of a; G R^, which 
seems essential since the Fourier transform was used in the definition. However, by 
means of local coordinates one may define pseudo-differential operators on a manifold 
(for example, the surface S), and then the principal symbol is independent of the choice 
of coordinates if ^ is interpreted covector. 

A basic (albeit non-obvious) property of pseudo-differential operators is that if 
P and Q are pseudo-differential operators of orders m, I, respectively, then their 
composition PQ is a pseudo-differential operator of order m + I, with principal symbol 
<Jm+i{PQ) = o'm{P)o'i{Q)- The identity operator has order with principal symbol 
cro(l) = 1, therefore if P is invertible then the principal symbol of its inverse is 

a_„(P-i) = MP)y'- 

In order to prove the shape- independent convergence of the high order surface modes 
we need one more well-known result (see [IE]): The Dirichlet-Neumann operators D± 
are pseudo-differential operators on S whose principal symbols are the same as in the 
case of a half-space (fT2l) and (|T3l) : 

ai(D_) = 1^1, ai(D+) = (16) 

From this we get a^i^D^) = — -j^ and then 

aoiD^D-') = -1. (17) 

This implies that the zeroth order operator R := D^D^^ + 1 has the principal symbol 
cro(-R) = —1 + 1 = 0, hence is in fact a pseudo-differential operator of order —1. Now 
any pseudo-differential operator of negative order on a bounded surface is a compact 
operator, and the spectral theory of compact operators implies that the eigenvalues of 
R form a a sequence converging to zero as k oo. The eigenvalues of D^D^ = —1+R 
are — 1 + r^, so the resonances obey 

efc = (—1 + rk)^^ ^ — 1 as A; — oo. (18) 

I Sometimes one considers even more general symbols, and then the ones we defined are called classical 
symbols of type (1,0), see [16]. 



Surface plasmon resonances of an arbitrarily shaped nanoparticle 



7 



We emphasise that this hmit is independent of the shape of the particle and equal 
to the result for the half space. The physical expectation based on our consideration of 
curvature becoming invisible to leading order at small scales is reflected in this argument 
by the fact that the principal symbol of the Dirichlet-Neumann operator is independent 
of the shape. 

We remark that with more refined methods of spectral asymptotics for pseudo- 
differential operators a more precise asymptotic result can be obtained which explains 
the correction in the formula = — 1 — for the sphere, ( |T4|) . 

Finally, we mention a different reduction of the plasmonic eigenvalue problem ([1])- 
([3]) to an eigenvalue problem on the surface S, derived by the method of layer potentials 
in |12], [13]: For r, r' on S let 



where n(r) is the outer unit normal. Then e = ^— ^ for the eigenvalues A in the equation 



Here a is the outward normal derivative (9„0 (which is proportional to the induced 
surface charge density). The operator on the left is compact (in fact, a pseudodifferential 
operator of order —1), hence 0, and this shows —1 again. For a proof of 

this fact see ^6] where also the relation of this operator and the Dirichlet-Neumann 
operators is discussed. 

5. Conclusion 

In this paper, we have applied methods of microlocal analysis to the study of plasmon 
resonances of arbitrarily shaped nanoparticles. We have first reformulated the boundary 
value problem for such resonances as an eigenvalue problem on the particle's surface. 
The fact that the Dirichlet-Neumann operators, which occur naturally in this context, 
are pseudo-differential operators then allows one to take advantage of the rich amount 
of knowledge available for these objects, and thereby to analyse the properties of the 
surface modes in rigorous mathematical terms. The remarkable ease with which the 
well-known result e = — 1 for the surface modes of a half-space has been recovered 
here is a clear indication for the power of this approach. Moreover, we have used the 
eigenvalue equation for proving that the limit of the high order modes is independent of 
the shape of the particle. As expected on the grounds of intuitive physical arguments, 
the high-order modes converge locally to the half-space modes. While this result itself 
appears natural, what matters here is the mathematical toolbox by which it has been 
obtained, which may be still somewhat unfamiliar to physicists, but which offers great 
conceptual clarity and flexibility. Thus, we hope that the approach suggested in this 
paper will prove useful for obtaining further insight into mathematical problems ocurring 
in plasmonics. 
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